A Harris- Todaro Agent-Based Model to Rural-Urban Migration 
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The Harris- Todaro model of the rural-urban migration process is revisited under an agent-based 
approach. The migration of the workers is interpreted as a process of social learning by imitation, 
formalized by a computational model. By simulating this model, we observe a transitional dynamics 
with continuous growth of the urban fraction of overall population toward an equilibrium. Such an 
equilibrium is characterized by stabilization of rural-urban expected wages differential (generalized 
Harris- Todaro equilibrium condition), urban concentration and urban unemployment. These classic 
results obtained originally by Harris and Todaro are emergent properties of our model. 



I. INTRODUCTION 

In this paper we turn upon the seminal Harris and To- 
daro Ijl work, which together with Todaro is consid- 
ered one of the starting points of the classic rural-urban 
migration theory 3]. The hypothesis and predictions of 
Harris- Todaro model have been subjected to econometric 
evaluation and have been corroborated by several studies 
0, 0, 0. • The key hypothesis of Harris and Todaro are 
that migrants react mainly to economic incentives, earn- 
ings differentials, and the probability of getting a job at 
the destination have influence on the migraton decision. 
In other words, these authors posit that rural-urban mi- 
gration will occur while the urban expected wage 
exceed the rural wage. From this crucial assumption, as 
denominated by Harris- Todaro 0, is deduced that the 
migratory dynamics leads the economic system toward 
an equilibrium with urban concentration and high urban 
unemployment . 

In our previous works 0, 0| we analyzed the rural- 
urban migration by means of an agent-based computa- 
tional model taking into account the influence of the 
neighborhood in the migration decision. The inclusion 
of the influence of neighbors was done via an Ising like 
model. The economic analogous to the external field in 
the Ising hamiltonian was the differential of expected 
wages between urban and rural sectors. Therefore, in 
theses works 8, 9| the crucial assumption of Harris and 
Todaro were taken for granted. 

Now, we are motivated by the following question: can 
the crucial assumption and equilibrium with urban con- 
centration and urban unemployment obtained from the 



original Harris- Todaro model be generated as emergent 
properties from the interaction among adaptative agents? 
In order to answer this question we implemented an 
agent-based computational model in which workers grope 
for best sectorial location over time in terms of earnings. 
The economic system simulated is characterized by the 
assumption originally made by Harris and Todaro. 

The paper is arranged as follows. Section ITTl describes 
the analytical Harris- Todaro model showing its basic 
equilibrium properties. In Section [llll we present the im- 
plementation of the computational model via an agent- 
based simulation and compare its aggregate regularities 
with the analytical equilibrium properties. Section IIVI 
shows concluding remarks. 



II. THE HARRIS-TODARO MODEL 

A. Assumptions 

Harris and Todaro 1] studied the migration of work- 
ers in a two-sector economic system, namely, rural sector 
and urban sector. The difference between these sectors 
are the type of goods produced, the technology of pro- 
duction and the process of wage determination. The ru- 
ral sector is specialized in the production of agricultural 
goods. The productive process of this sector can be de- 
scribed by a Cobb-Douglas production function: 



(1) 
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where Ya is the production level of the agricultural good, 
Na is the amount of workers used in the agricultural pro- 
duction, Aa > and < < 1 are parametric constants. 

Similarly, the urban sector also has its productive pro- 
cess described as Cobb-Douglas production function: 



(2) 



where Ym is the production level of the manufactured 
good, N„i is the quantity of workers employed in the 
production of manufactured goods, Am > and < a < 
1 are parametric constants, jlij 

Both goods and labor markets are perfectly compet- 
itive. Nevertheless, there is segmentation in the labor 
market due to a high minimum urban wage politically 
determined. 

In the rural sector, the real wage, perfectly flexible, is 
equal to the marginal productivity of labor in this sector 



(3) 



where Wa is the real wage and p is the price of the agri- 
cultural good, both expressed in units of manufactured 
good. 

In the urban sector, a minimum wage, Wm, is assumed 
fixed institutionally at a level above equilibrium in this 
labor market. It can be formalized as 



such that Nm < Nu, 



(4) 



where Nu is the amount of workers in the urban sector. 

The relative price of the agricultural good in terms of 
the manufactured good, p, varies according to the relative 
scarcity between agricultural and manufacturated goods. 
Then, 
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B. Temporary Equilibrium 

Given a parametric constant vector 
(Aq, ^,„, 0, a, p, 7), an initial urban population iV„, 
and a minimum wage Wm one can calculate the tempo- 
rary equilibrium of the economic system by using eqs. 
(ITlHt . 

From eq. I@J one can find the employment level at the 
manufacturing sector 



N„ 



f aAr, 



(7) 



Replacing eq. Q in eq. (O we get the production level 
of the manufacturing sector 



Yra = A^r° . (8) 

From eq. © one can obtain the relation 



Nn=N- N.,, 



(9) 



which is used with eq. 1^ to obtain the agricultural 
production 



Ya =Aa{N- NuY 



(10) 



By using eqs. ijHJl and H1U|) the terms of trade are 
determined 



P = P 



"^"^ ° (^) 
Aa{N - Nu)^ 



(11) 



Finally, by using eqs. (O, JHl and the rural wage 
in units of manufacturated good is obtained 



where p > Q and 7 > are a parametric constants. 7 is 
the elasticity of p with respect to the ratio l^/F^. 

The overall population of workers in the economy is 
iV, which is kept constant during the whole period of 
analysis. By assumption there are only two sectors and 
rural prices are wholly flexible, which implies that there 
is full employment in the rural area, i.e., all workers living 
at the rural sector are employed at any period. Then at 
any period the following equality is verified: 



Na+N^ = N. (6) 



In sum, the vector {N„i, Ym, Na, Ya,p, Wa) configures a 
temporary equilibrium that might be altered whether oc- 
curs a migration of workers, induced by the differential of 
sectorial wages, which changes the sectorial distribution 
of overall population. 

C. The Long Run Equilibrium 

Harris and Todaro, in determining the long run equi- 
librium, i.e., the absence of a net rural-urban migratory 
fiow, argue that the rural workers, in their decision on 
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migrating to the urban area, estimate the expected ur- 
ban wage, w!:, defined as: 



Nrr, 



(13) 



The ratio Nm/Nu, which is the employment rate, is an 
estimative of the probabihty that a worker Hving at urban 
sector gets a job in this sector. 

As mentioned before, the key assumption of the model 
of Harris and Todaro is that there will be a migratory flow 
from the rural to the urban sector while the expected 
urban wage is higher than the rural wage. Thus, the 
long run equilibrium is attained when the urban worker 
population reaches a level such that the expected urban 
wage equates the rural wage: 



Wa 



0. 



(14) 



This equality is known in the economic literature as 
the Harris- Todaro condition. Harris and Todaro argue 
that the differential of expected wages in eq. H14|l can be 
a constant value (5^0. When this differential reaches S, 
the net migration ceases. This generalized Harris-Todaro 
condition can be expressed as follows: 



-Wa= 5. (15) 

The level of the urban population that satisfies the eq. 
((T3)l . i.e., the equihbrium urban share n* = N*/N, is 
determined from the solution of the equation resulting 
from substitution of equations H12I) . H13|) in eq. (|15|l : 



Nu \Wm 



{N - N^y-'t'+'t't 



^5. 



(16) 



The solution of eq. (|16|l is parametrized by the vector 

m: P: 0i ^ni)- 

Harris and Todaro ^j] , in order to evaluate the stability 
of the long run equilibrium, postulate a mechanism of 
adjustment that is based on the following function of sign 
preservation: 



iV„ = V«-Wa), with V'(-) > and ipiO) = 0. (17) 

The differential equation that governs the state tran- 
sition in the model of Harris and Todaro is obtained by 
replacing equations ((T^ . ifT^ in eq. ifTTji . Based on 
this postulated adjustment process, Harris and Todaro 
Q show that the long run equilibrium is globally asymp- 
totically stable. This means that the economy would tend 



FIG. 1: Numerical solution of eq. for different values of 
Wm- Squares; urban share n* ; Circles; urban unemployment 



rate (1 — Nm/N). Fixed parameters used are Aa 
1.0, (j) = 0.3, a = 0.7, p = 1.0 and 7 = 1.0. 
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FIG. 2; Numerical solution of eq. Hl()|l for different values of 7. 
Squares; urban share nJl; Circles; urban unemployment rate, 
(1 — N*/N). Fixed parameters used are Aa — 1.0, Am = 1.0, 
(j> = 0.3, a = 0.7, p = 1.0 and Wm = 0.8 



to long run equilibrium with unemployment in the urban 
sector generated by the presence of a relatively high min- 
imum wage for all possible initial conditions. From now 
on we will refer to the long run equilibrium simply as 
equilibrium. 

Based on the numerical solutions of eq. (|16() one can 
evaluate the impact that the variation of the minimum 
wage and the elasticity of the terms of trade on the equi- 
librium. In Figure ^ we see that under the hypothesis of 
a Cobb-Douglas technology, the equilibrium urban share, 
n*, does not depend on the minimum wage Wm- How- 
ever, changes in the value of Wm reduces the labor de- 
mand on the manufacturing sector what results in higher 
unemployment rates in the equilibrium. 

In turn, as seen in Figure |21 changes in the elasticity 
of the terms of trade alter slightly the equilibrium urban 
share and unemployment rate. A net migration toward 
urban sector shift the terms of trade to higher values. 
The greater 7 the greater this shift, what cause an in- 
crease in the rural wage in units of manufacturing good. 
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becoming the urban sector less attractive. 



III. HARRIS-TODARO AGENT-BASED MODEL 

In this section we describe the implementation of the 
computational model we proposed, as well as the aggre- 
gate patterns obtained numerically and the comparison 
with the respective analytical results. 



A. Computational Implementation 

Initially, workers are randomly placed in a square lat- 
tice with linear dimension L = 500. The reference values 
of the parameters used for these simulations are the same 
done to evaluate the equilibrium of the Harris- Todaro 
model, namely, Aa = 1.0, A,n = 1.0, = 0.3, a = 0.7, 
p = 1.0 and 7 = 1.0. The value of the minimum wage 
used is Wm = 0.8 and the initial urban fraction of the to- 
tal population is n„ = 0.2, where n„ = N^/N is the nor- 
malized urban population also called urban share. The 
initial value riu = 0.2 is in agreement with historical data 
of developing economies. Given these parameters, one 
can calculate the vector which characterizes temporary 
equilibrium of the system by using eqs. (|7I12I) . 

By using eq. (0, the employment level of the urban 
sector, Nrtn is obtained. If n„ < Nm/N all workers in the 
urban sector are employed and each individual i earns the 
wage given by the manufacturing marginal productivity, 
Wi = aAmN"^^ . Otherwise, Uu > Nm/N there will be 
a fraction of Nm/Nu workers employed, which earn the 
minimum wage, Wi = Wm, and (1 — Nm/Nu) workers 
unemployed, which earn a wage Wi = 0. 

Each worker can be selected to review his sectorial lo- 
cation with probability a, called activity (Tlj. Therefore, 
in each time step only a fraction of workers becomes po- 
tential migrants, going through the sectorial location re- 
viewing process. Potential migrants will determine their 
satisfaction level of being in the current sector by com- 
paring their earnings, Wi, among nearest neighbors. 

The potential migrant starts the comparison process 
with a initial satisfaction level Si — 0. When Wi > 
Wneighbor the Satisfaction level st is added in one unit; 
if Wi < Wneighbor, Si IS diminished in one unit; if Wi ~ 
Wneighbor, Si docs uot changc. After the worker has 
passed through the reviewing process his/her satisfaction 
level is checked. The migration will occur only if Si < 0, 
what means that the worker's i earnings is less than the 
most of his/her nearest neighbors. 

After all the potential migrants complete the reviewing 
process and have decided migrate or not, a new config- 
uration of the system is set. Therefore, once again a 
new temporary equilibrium of the system is calculated 
by using eqs. (|8I12|I . The whole procedure is repeated 
until a pre-set number of steps is reached. It is impor- 
tant to emphasize that Nm is kept constant throughout 
the simulation. Its given by eq. Q which depends on 



the technological parameters, a, Am, and the minimum 
wage, Wm, which are constants too. 



B. Analysis of the Emergent Properties 

In this section we develop the analysis of the long run 
aggregate regularities of Harris- Todaro agent-based com- 
putational model. These long run properties will be com- 
pared between the solution of the analytical model and 
simulations we ran. 




FIG. 3: Urban share function of simulation steps. From 

top to bottom the initial urban shares are 0.9, 0.8, 0.3, 0.2. 

Figures Ol El and O show the basic characteristics of the 
transitional dynamics and long run equilibrium generated 
by simulations. When the economic system has a low 
initial urban share, n„ = 0.2 or n„ — 0.3, there is a 
net migration toward urban sector. This migration takes 
the urban sector from a full employment situation to an 
unemployment one. The positive differential of expected 
wages that pulls workers to the urban sector diminishes. 
However, if the economic system initiates with a high 
urban share, Uu — 0.8, or n„ — 0.9 there is net flow 
of migration toward rural sector in such a way that the 
unemployment rate of the urban sector decreases. In this 
case, the differential of expected wages is negative. 

In an economy mainly rural (ri„ < 0.5), the transi- 
tional dynamics characterized by a continuous growth 
of population of the urban sector with a differential of 
expected wages relatively high is followed by the stabi- 
lization of rural-urban differential of expected wages. In 
other words, the generalized Harris- Todaro condition, eq. 
(|15|l . arises as a long run equilibrium result of the agent- 
based migratory dynamics. 

FigureOlalso shows that even after the urban share has 
reached an stable average value, there are small fluctua- 
tions around this average. Therefore, differently from the 
original Harris- Todaro model, our computational model 
shows in the long run equilibrium the reverse migration. 
This phenomenon has been observed in several develop- 
ing countries as remarked in Ref. 

In Figures El [3 and |S1 one can see that for a given value 
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FIG. 4: Unemployment rate (1 — Nm/Nu) as function of sim- 
ulation steps. From top to bottom the initial urban shares 
are 0.9, 0.8, 0.3, 0.2. 




FIG. 6: Equilibrium urban share riu as function of the tech- 
nological parameter a and the minimum wage Wm- White 
area is not a valid combination of parameters. 
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FIG. 5: Rural-urban expected wage differential (w^ — Wa) as 
function of simulation steps. From top to bottom the initial 
urban shares are 0.9, 0.8, 0.3, 0.2. 



of a, the variation of Wm practically does not change the 
equilibrium values of the urban share, the differential of 
expected wages and the unemployment rate. However, 
for a given Wm, higher values of a make the urban sector 
less attractive due the reduction of the employment level. 
This causes a lower equilibrium urban share, a higher 
unemployment rate and a gap in the convergence of the 
expected wages. 

In Figures El El ^^^d ^] can be seen that for a fixed 
value of 7, the equilibrium values of the urban share, the 
differential of expected wages and unemployment rate do 
not have a strong dependence with Wm- However, varia- 
tions in 7 for a fixed Wm, dramatically change the equi- 
librium values of the variable mentioned before. Higher 
values of 7 generate a lower urban concentration, a higher 
gap in the expected wages and a higher unemployment 
rate in the equilibrium. 

Finally, in Figure [T^ is shown that the convergence of 
migratory dynamics for a urban share, compatible with 
historical data, is robust in relation to the variation of 
the key technological parameters, a and (j). 



FIG. 7: Equilibrium differential of expected wages as function 
of the technological parameter a and the minimum wage Wm ■ 
White area is not a valid combination of parameters. 
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FIG. 8: Equilibrium urban unemployment rate (1 — Nm/Nu) 
as function of the technological parameter a and the min- 
imum wage Wm- White area is not a valid combination of 
parameters. 
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FIG. 9: Equilibrium urban share Tifi as function of the pa- 
rameter 7 and the minimum wage Wm- White area is not a 
vaUd combination of parameters. 



FIG. 12: Equihbrium urban share Tiu as function of the tech- 
nological parameters a and 4>- White area is not a vahd com- 
bination of parameters. 
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FIG. 10: Equihbrium differential of expected wages (wZ — Wa) 
as function of the parameter 7 and the minimum wage Wm- 
White area is not a valid combination of parameters. 



FIG. 13: Equilibrium differential of expected wages (wZ — Wa) 
as function as function of the technological parameters a and 
(p. White area is not a valid combination of parameters. 
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FIG. 11: Equilibrium urban unemployment rate {1 — Nm/Nu) 
as function of the parameter 7 and the minimum wage Wm- 
White area is not a valid combination of parameters. 



FIG. 14: Equilibrium urban unemployment rate {l~Nm/Nu) 
as function of the technological parameters a and 4>. White 
area is not a valid combination of parameters. 
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IV. CONCLUSION 

In this paper we developed and agent-based computa- 
tional model which formalizes the rural-urban allocation 
of workers as a process of social learning by imitation. 
We analyze a two-sectorial economy composed by adap- 
tative agents, i.e., individuals that grope over time for 
best sectorial location in terms of earnings. This search 
is a process of imitation of successful neighbor agents. 

The dispersed and non-coordinated individual migra- 
tion decisions, made based on local information, generate 
aggregate regularities. Firstly, the crucial assumption of 
Harris and Todaro, the principle that rural-urban migra- 
tion will occur while the urban expected wage exceed the 
rural wage, comes out as spontaneous upshot of interac- 
tion among adaptative agents. 

Secondly, the migratory dynamics generated by agents 
that seek to adaptate to the economic environment that 
they co-create leads the economy toward a long run equi- 
librium characterized by urban concentration with ur- 
ban unemployment. When this long run equilibrium is 



reached, the generalized Harris-Todaro condition is sat- 
isfied, i.e., there is a stabilization of the rural- urban ex- 
pected wage differential. 

Thirdly, the impact of the minimum wage and elastic- 
ity of terms of trade in a long run equilibrium obtained by 
simulations are in agreement with the predictions of the 
original Harris-Todaro model with Cobb-Douglas tech- 
nology. 

Finally, the simulations showed an aggregated pattern 
not found in the original Harris-Todaro model. There is 
the possibility of small fluctuations of the urban share 
around an average value. This phenomenon is known as 
reverse migration. 
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